
Rapid Neutron Capture in Supernova Explosions 


I 


/ 


J.W. Truran* and W.D. Arnett** 

— Institute for Space Studies 
Goddard S p a gg F L i -ght- Ce nteje.,- NASA 
New York, N.Y._i 



S . Tsuruta 

Smithsonian Astrophysical Observatory 
Cambridge, Mass. 


and 


A.G.W. Cameron / 

Belfer Graduate School of Science 
Yeshiva University 
New York, N.Y. 


and 

Institute for Space Studies 
Goddard Space Flight Center, NASA 
New York, N.Y. 



§ 

\ 

\ 










o 

U. 



co 

X 

5 



UJ 

> 

© 


UJ 

o 

a. 



y 

cc 

o 

a 

O 

in 

S 

CL 

TJ 

o 

Urn 

cc 

0. 

o 

H 

if) 

w 

CO 

X 

u 

i 

3 

-> 

n 

1.0 

Q. 

LL 



© 

o 

o 



!fc 


Abstract 


The implications of recent studies of the dynamics of the cores 


of highly evolved massive stars are considered with regard to the 
general problems of nucleosynthesis. The typical conditions estimated 
for these models are shown to be very promising for the process of 
element synthesis by neutron capture on a fast time scale (the r- 
process of Burbidge et al, 1957 )- 
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I. Introduction 

Studies of the thermonuclear reaction mechanisms believed to be 
responsible for the production of the heavy elements (A > 40) have 
generally been restricted by the lack of a realistic model of the 
stellar environment. Investigations of both the slow and rapid neutron 
capture processes have involved attempts to fit the observed abundance 
patterns with various choices for the available neutron flux (Seeger 
et al , 1965 )* Considerations of the formation of nuclei in the iron 
peak region have generally been concerned with the deviations of the 
observed abundances flora an equilibrium configuration (Burbidge, 

Burbidge, Fowler and Hoyle, 1957) or with the nuclear transformations 
which take place under somewhat arbitrarily specified temperature and 
density conditions (Truran et al, 1966). While such studies have 
established the role of these processes in the synthesis of the heavy 
elements, they have not provided a complete picture of the sequence of 

these processes within the general framework of stellar evolution and 

» 

nucleosynthesis . 

* Recent studies of the dynamics of the cores of highly-evolved massive star 
(Arnett, 1966 and 1967a; Colgate and White, 1966) have revealed that for cores 
of mass M c ^ 4 a considerable fraction of the core mass will be ejected. 
This mass ejection may be identified with supernova explosions. At 
the extreme temperature and density conditions realized in the core the 
matter is composed predominantly of neutrons. As the ejected matter 
expands and cools, these neutrons may be converted to protons, helium 


and heavier elements. Preliminary estimates indicate that element 
synthesis may take place readily under these conditions. In this paper 
we will consider the implications of this model with regard to the 
problem of the synthesis of the heavy elements by rapid neutron capture. 
The possible production of heavy elements in the supernova envelope is 
also discussed. The formation of iron peak elements in the supernova 
envelope has been considered in another paper (Truran et ^1, 1967). 

The outline of this paper is as follows. In Section II we review 
in detail the supernova model on which our calculations are based. In 
Section III the nuclear reaction network employed in these calculations 
is defined and a brief discussion of the relevant nuclear physics is 
presented. The results of our studies of the extent of heavy element 
production are presented in Section IV. In Section V the abundance 
configurations predicted in equilibrium for our temperature and density 
conditions are presented. The three processes of element synthesis 
which are believed to be responsible for the production of the heavy 
elements are defined in Section VI. In Section VIII, the implications 
of the supernova conditions for the process of element synthesis by rapid 
neutron capture are considered in detail. A brief discussion of other 
possible neutron capture processes is given in Section VIII. 

II. Supernova Model 

The paths in the density-temperature plane of two representative 
mass zones of the 2 supernova core model of Arnett (1967a; see also 
Arnett and Cameron, 1967) are shown in Figure 1. In this paper we will 
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be concerned with the innermost ejected layers of the core, as represented 
by the mass zone at O .9743 • The initial model evolves along the 

adiabat AB because hydrodynamic compression has become rapid compared to 
neutrino energy loss. The evolution of the 0.9743 ^ interior mass 
zone is extremely rapid from B to C. Inadequate sampling of the calcula- 
tion (once per thousand time cycles) makes our knowledge of the details 
of the path uncertain here. Evolving from C to D this material encounters 
the thermal wave due to energy-transporting neutrinos. After heating, 
the infall of the material is reversed and it begins to expand from 
D to E. At E it again encounters the neutrino emission surface. The 
nonadiabatic evolution along EF is in the thin- approximation region for 
neutrino transport. The heating processes exceed cooling because the 
core is expanding and therefore decreasing its "optical" depth for 
neutrinos (antineutrinos); this results in more rapid energy transfer 
from the core. Along FG the cooling is almost adiabatic because the 
neutrino flux is attenuated by spherical geometry and the evolution is 
too fast for the usual pair annihilation neutrino loss rate to be 
important (although it is included in the calculation) . At G the equi- 

4 

librium for strong interactions transforming He to nucleons begins to 

4 

favor He again. 

In the calculation shown in Figure 1, the ratio of neutron and 
proton number densities was of order unity, so that evolution proceeded 
along GH. This was an oversimplification of the time-dependent equation 
of state. Electron capture was assumed to occur only on nuclei, and 
only if the density was more than 10“'“'*' gm.cm. 


While this assumption 
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does not affect the behavior of the material which forms the nuclear- 
density central region, and has little effect on the qualitative hydro- 
dynamic behavior of the ejected mass, a more accurate estimate of the 
neutron to proton ratio is desirable for a discussion of the implications 
of this model for element synthesis . 

In a subsequent consideration of this problem (Arnett and Cameron, 

196 7) it was found that the weak interactions should freeze out in the 

early stages of expansion of the ejected material at a temperature of 

9 -3 

approximately 20 billion degrees and a density of 7x10 gm.cm. : the 
neutron to proton ratio realized in this freezing process was approxi- 
mately 8. This ratio was not found to be a sensitive function of the 
freezing density. For this neutron to proton ratio only a fraction of 
the mass will be converted to helium; energy considerations suggest that 
the evolutionary path shown in Figure 1 will be amended in the manner 
Indicated by the dashed line . 

The temperature -density history of the stellar material employed In 
our considerations of nucleosynthesis in this paper are those obtained 
in the manner described above. While these results are the best estimates 
available to date, it should be emphasized that a proper coupling of 
the hydrodynamic equations, the transport equations for weak interactions 
and compositional changes has not yet been accomplished, and that 
rotational effects have not been included (see Colgate, 1967 and Arnett, 
1967b) . Under these circumstances it seems most prudent to consider 
that a range of adiabats may be involved either in a single supernova 
or among a group of supernovae . 


5 


The temperature -density profiles employed in our calculations are 

shown in Figure 2. The initial temperature we consider is lO 1 ^ °K, at 

which point the weak interactions have frozen out and the formation of 

helium from neutrons and protons will proceed rapidly. The density 

curve shown here is that which has been corrected for the neutron rich 

behavior; this is the low adiabat . For the high adiabat, the density is 

reduced by a factor of 10. The decrease in temperature and density with 

time results from the expansion and cooling of the ejected material. 

We are interested in the subsequent thermonuclear evolution of the 

matter ejected on these supernova models . For the temperature and 

densities of interest here, the time scale for the conversion of neutrons 
4 

and protons to He is considerably less than the hydrodynamic time scale. 
The critical consideration is whether the buildup toward heavier nuclei 
will continue, providing seed nuclei for the process of neutron capture 
on a fast time scale . 

III. Nuclear Reaction Network and Reaction Rates 

It was shown in a previous paper (Arnett and Cameron 1967 ) that 
12 

substantial amounts of C are built up in equilibrium with neutrons, 
protons, and alpha-particles during the expansion of neutron-rich material 
along the adiabats discussed above. We wish to determine whether a true 
nuclear statistical equilibrium will be approached and maintained during 
the expansion. The nuclear reaction pathways along which the approach 
to equilibrium will take place presumably lie in the region of very 
neutron-rich heavy nuclei. The approach to nuclear statistical equilibrium 
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has so far been investigated only for pathways which lie close to the 
ordinary valley of beta stability (Truran, Cameron, and Gilbert 1966). 

We do not feel that the techniques used to determine the appropriate 
nuclear reaction cross sections (Truran, Hansen, Cameron, and Gilbert 
1966) will give reliable predictions for very neutron-rich nuclei, and 
hence extensive calculations have not been made with such nuclei. 

Instead, we have adopted the point of view that if nuclear statistical 
equilibrium can be approached by any set of reaction pathways during 
the expansion along the adiabat, then it should certainly be attained by 
the fastest pathway. We have been conservative in that we have required 
that charged-particle addition reactions must involve alpha-particle 
capture at least up to nuclear charge numbers equivalent to those encountered 
in nuclear statistical equilibrium. We have chosen a specific pathway 
lying close to the ordinary valley of beta stability so that the reaction 
rates should be fairly reliable. 

Accordingly, we have defined a simple nuclear reaction network 

providing suitable nuclear reaction links connecting neutrons, protons 

86 

and alpha-particles to Kr. This network is illustrated in Figure 3 * 

In the light element region we have included the more important reactions 
linking neutrons, protons and alpha -particles* as indicated in the figure. 

The rates employed for these reactions were taken from Wagoner, Fowler 
and Hoyle (1967). 

For the heavy elements, we have considered only alpha-particle and 
neutron capture reactions. Where experimental date is available, these 
rates are determined as sums over individual resonance contributions; 
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otherwise, they are calculated theoretically from average nuclear 
properties (Truran et al, 1966a). Specifically, experimental determi- 
nations of the reaction rates were employed for the triple-alpha reaction 

12 16 20 24 

and for the alpha-capture reactions on C, 0 , Ne, and Mg 
(Truran et al, 1966a, b; Reeves 1965 ; Fowler et al 1967). 

Shaw and Clayton (1967) have considered the enhancement of radiative 
transition rates by inelastic scattering with electrons and ions in 
stellar matter. They have pointed out that these effects may be important 
under non-equilibrium conditions in fusion reactions in the cores of 
supernovae. In our calculations, both carbon and the heavy elements 
remain in equilibrium with respect to neutrons and alpha-particles down 
to a temperature of approximately five billion degrees. At this tempera- 
ture, the density for the low adiabat considered in our calculations is 
7 -3 

8 . 75 x 10 gm.cm. and clearly those effects are not important. In 
general, this process will become important in mass zones characterized 
both by considerably higher densities and by shorter expansion time 
scales . 

The relaxation of this nuclear reaction network is governed by the 
equations for the rates of change of the abundances of the various 
nuclear species. The non-linear equations for the changes in the number 
densities of the various constituents have been expressed in a backward 
difference form and solved by matrix inversion. 
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IV. Production of Seed Nuclei 


We have considered the nuclear transmutations which will take place 

for a range of temperature and density conditions. In Figure 4 are 

plotted the mass fractions of the more important nuclear constituents 

as a function of time for the high density profile. The initial neutron 

to proton ratio is 8. It is evident from the behavior in the early 

stages that an equilibrium has been rapidly established among neutrons, 

12 

protons, alpha-particles and C. As the temperature decreases, and 

12 

the alpha -particle capture past C becomes more rapid, the abundances 

12 

of the light elements and C decreases while the abundances of inter- 

86 

mediate nuclei and ultimately Kr are rapidly increasing. At a 

-2 

temperature of approximately 4 billion degrees (approximately 2x10 sec.) 

the reactions involving the heavy elements freeze out in the sense that 

the nuclear reaction time scale becomes long compared to the expansion 

86 

time scale. The final abundance of °°Kr is 20$ by mass. 

The results of a similar calculation performed for the lower density 

12 

profile are shown in Figure 5. Here the build-up of C is somewhat 

86 

impeded by the lower density. However, the final Kr abundance still 
approaches 1 QPjo by mass . 

The effects of varying the total neutron to proton ratio realized 

in freezing are shown in Figure 6. Here the ratio N/P=3, and the low 

density adiabat is employed. Generally, the protons will be converted 

completely to helium, hence a larger fraction of the mass can be 

86 

converted to heavy elements. The Kr abundance realized in this calcu- 
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lation approaches 50 % by mass. 

For both of the density profiles considered above, an appreciable 

86 

fraction of the mass was found to be processed to Kr. It is of 
interest to determine more precisely the limiting conditions under which 
this buildup process will take place . In Figure 7 are shown the results 
of a similar calculation performed for an order of magnitude higher 
adiabat : 

P = 7xl0 3 T 3 

7 

(This is a somewhat unrealistic condition with regard to the supernova 

model discussed previously; further, the neutron to proton ratio employed 

is rather high for these conditions.) The abundance of heavy elements 

46 

produced is approximately five percent by mass . Here Ca is the 
most abundant nucleus, as heavy element production is impeded by the 
low densities. 

12 

These results suggest that the breakthrough past C to heavy 
elements can take place readily under these conditions. It is clear, 
however, that our simple linear chain past carbon does not provide an 
accurate representation of the flows through the iron peak and beyond. 
Specifically, for the high neutron concentration in our calculation the 
major flow might proceed far into the neutron rich regions. In this 
regard, these calculations should give an upper limit both on the time 
scale required and on the ratio of neutrons to seed nuclei for rapid 


neutron capture . 
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V. Equilibrium Abundances 


The final distribution of seed nuclei for the rapid neutron capture 
process may be estimated on the assumption that the nuclear reactions 
are sufficiently rapid that a condition of nuclear statistical equilibrium 
is established (Hoyle, 1946; Tsuruta and Cameron, 1965). The abundance 
of a nucleus of atomic number Z and mas number A in equilibrium is 
given by 


n(A,Z) = w(A,Z) 


/ 2nh 2 \ 
\MkT J 


| (A-l) 


n (A-Z) n. Z 

(A 3 / 2 ) — exp[Q(A,Z)/kT] 

^A 


where n and n are the number densities of free protons and neutrons 
p n r 

and Q(A,Z) and w(A,Z) are the binding energy and partition function, 
respectively, of the nucleus (A,Z). 

These equations are supplemented, for our conditions, by the equation 
for the conservation of mass 


P = l A k W/K 0 

and by an equation specifying the total ratio of neutrons to protons 

“A - 1 <VV n <\’ z k>/! W 

Here N o is Avagadro's number, and the sums are taken over all nuclear 
constituents including neutrons and protons. If we specify the tempera- 
ture, the density and the tota l ratio of neutrons to protons these equations 
can be readily solved for the equilibrium number densities. 
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The equilibrium abundances obtained for a temperature of 4x10^ °K, 
a density of 4.48x10 gm.cm. ~ > and a total neutron to proton ratio of 
eight are shown in Figure 8. These conditions are those predicted by 
the low adiabat near the point at which the reactions involving the 
heavy elements have frozen out. The masses employed for nuclei in the 
neutron rich regions far from the valley of beta stability are those 
predicted by the exponential mass formula of Cameron and Elkin ( 1965 ). 

iyO 

The peak nucleus in equilibrium under these conditions is Ni, 

having 28 protons and 50 neutrons. It is understandable that this 

nucleus, possessing "magic" numbers of both protons and neutrons, is 

exceptionally stable. However, at slightly lower temperatures the peak 
78 120 

shifts from Hi to Sr (38 protons and 82 neutrons). This behavior 

is evident in Figure 9j where the equilibrium abundance distribution 

9 o 

for the same adiabat at a temperature 3x10 K is shown. The major 
abundance peak under these conditions is thus seen to be one neutron 
closed shell removed from the valley of beta stability. 

The equilibrium abundances for the high adiabat at the same freezing 
temperature are shown in Figure 10. The total ratio of neutrons to 
protons is again taken to be 8. At the factor of ten lower density 

rv Q 

considered here, the abundance peak remains in the vicinity of Ni for 

Q Q 

a temperature of 3x10^ K. The shifting of the peak to the N = 82 

9 o 

neutron shell occurs for T — 2xKr K, well below the temperatures for 
which any appreciable buildup of heavy elements by charged particle 
reactions can take place. 
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The temperature dependences of the abundances of neutrons, protons 

and alpha particles are shown in Figure 11. The total number density 

of nuclei with Z > 2 (denoted by 2 n(A,Z)) is also included. The 

Z>2 

solid lines correspond to the high density adiabatic cooling curve and 
the dashed lines to the low density case. The peak nuclei under 
different temperature -density conditions are indicated by their respective 
mass numbers and charge numbers (A^,Z^). 

The general trends evident in these equilibrium studies may be 
summarized as follows. A decrease in temperature , for the same adiabatic 
density dependence, results in a shifting of the abundance peak toward 
more neutron rich nuclei. The presence of increased abundances of 
neutron rich nuclei in equilibrium is also favored by higher densities. 
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VI. Mechanisms for Heavy Element Production 

There are three distinct processes of nucleosynthesis which are 
thought to be responsible for the production of the heavy elements 
(Burbidge et al, 1957; Cameron, 1957): neutron capture on slow and fast 
time scales and the process responsible for the production of the by- 
passed nuclei. The abundance trends of even mass numbers for these 
processes are shown in Figure 12 (Cameron, 1967). 

Neutron capture on a slow time scale (the s -process) is characterized 
by neutron capture lifetimes which are long compared to typical beta 
decay lifetimes in the vicinity of the valley of beta stability; the path 
of this process therefore proceeds along the valley of beta stability, 
the abundance distribution peaking at the neutron closed shell positions 
where the neutron capture cross sections decrease. These are evident 
in Figure 12 near mass number A = 90 ? 138 and 208. 

If the neutron capture lifetimes are short compared to the appropri- 
ate beta decay lifetimes, the neutron capture path will lie far out on 
the neutron rich side of the valley of beta stability. Such a situation 
might be realized in the presence of an extreme neutron flux; the tempera- 
ture must be sufficiently low that progress along the neutron capture 
chains is not impeded by photodisintegrations. Following the exhaustion 
of the neutron flux, the capture products approach the valley of beta 
stability by a series of beta decays. This sequence corresponds to 
neutron capture on a fast time scale (the r-process). 
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The abundance features near mass numbers A — 130 and 195 are 
generally attributed to this rapid neutron capture process (Seeger et al , 

1965 ). 

The bypassed nuclei are nuclei on the neutron deficient side of the 
valley of beta stability which cannot be produced by any neutron capture 
process. In principle , these can result from neutron photodisintegration 
processes or proton capture processes (the p-process) proceeding on 
the products of earlier neutron capture synthesis. The apparent 
constancy of the ratio of the abundance level of bypassed nuclei to 
that for the neutron capture products , evident in Figure 12, is consistent 
with this picture . 

In the subsequent discussion we will be concerned primarily with 
the conditions under which the production of heavy elements by rapid 
neutron capture can take place . 

VII . Rapid Neutron Capture 

We are now prepared to consider the implications of the supernova 

shock conditions with regard to the problem of element synthesis by 

rapid neutron capture. We have found that the production of elements 

in the vicinity of the iron peak and beyond can take place in the early 

stages of expansion of the ejected material via charged-particle and 

neutron capture reactions. The charged particle reaction rates decrease 

rapidly with temperature and a stable abundance configuration is realized 

* 9 o 

at a temperature in the vicinity of 4x10 K. This abundance pattern 
has been distorted severely by the limited nuclear reaction network we 
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employed; the abundances predicted in equilibrium for the "freezing” 
conditions should provide a more realistic estimate of those realized 
in the expanding material* 

While the charged-particle reactions are no longer contributing to 
the buildup of heavy elements below 4x10^ °K, and progress toward elements 
of higher Z is impeded , neutron capture buildup under these conditions 
should be limited only by neutron photodisintegration reactions. The 
equilibrium calculation displayed in Figure 10 predicts a neutron number 
density of 1.88xlO J cm. and a ratio of neutrons to heavy nuclei of 
approximately 200. As the temperature is further decreased , the neutron 
photodisintegration rates will decrease, allowing successive neutron 
captures to proceed toward nuclei which are increasingly beta-unstable. 

Seeger et_ al ( 1965 ) have pointed out that if the abundance peaks 
at A = 130 and 195 are to be understood in terms of nuclear shell 
structure, an environment characterized by a high temperature and neutron 
number density (T =* 10^ °K, n n 10^ cm. ~ > ) is required. These 
conditions are necessary, on their model, for the establishment of a 
thermal equilibrium between the isotopes of any element, neutrons and 
photons on a time scale which is short compared to typical beta-decay 
lifetimes. It is clear, however, that these conditions are extremely 
sensitive to mass formula predictions of the neutron binding energies 
and beta-decay energies of the neutron rich isotopes formed in the 
rapid neutron capture process . 

We feel that the rapid neutron capture buildup will be terminated 
either by the exhaustion of the neutron flux or by the expansion of the 
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medium to extremely low densities. In order to determine which of these 
is the limiting consideration, and to calculate the final abundances 
resulting from this neutron capture process, the subsequent history of 
the medium must be followed in detail. This will involve coupling the 
hydrodynamic equations with the equations governing the nuclear trans- 
formations. Two of us (W.D. Arnett and J.W. Truran) are currently 
involved in the formulation of this problem. 

Taken at face value, the conditions we have considered here lead 
to the formation of a large amount of heavy elements. It also seems 
likely that many of the elements present in nature in the silicon to 
iron region have been made in supernova explosions (Truran et al, 1967). 

The products of the rapid neutron capture process have an observed abundance 

-4 -6 

relative to the silicon to iron region in the range 10 to 10 . It 

is clear, therefore, that of all the mass processed in supernova 
explosions only a small part can be processed in the manner outlined 
in this paper. The trends in current supernova models are consistent 
with this conclusion. 

VIII . Other Neutron Capture Processes 

There are several features of the observed abundance distribution 
of the heavy elements which are difficult to account for on the basis of 
the rapid neutron capture process outlined above; specifically the 
abundance peak in the vicinity of mass number A ~ 103 and that in the 
rare earth region (see Figure 12 ). We believe these may have resulted 
from the exposure of a region containing a few percent by mass of 
the products of slow neutron capture to a somewhat smaller 
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neutron flux. This condition may be realized in the supernova envelope 

in the vicinity of the helium burning shell source, during the passage 

of the shock wave (Truran et al, 1967)- We have found that neutrons 

25 26 29 

released in (a ? n) reactions proceeding on Mg, Mg and y Si resulted 

l8 8 

in a neutron number density exceeding 2x10 particles/cm. on a time 

-2 

scale of the order of 10 seconds. In this case the maximum shock 

9 o 

temperature must lie close to 3x10 K, sufficient to produce rapid 
(o ,n) reactions but not sufficient to destroy the heavy nuclei by 
photodisintegration. Detailed calculations of the buildup of the 
anomalous fast capture peaks resulting from neutron capture on nuclei 
in the preceding slow neutron capture peaks near A = 90 and 140 have 
been found to be extremely sensitive to our estimates of the neutron 
capture cross sections. We are currently undertaking a reevaluation 
of these cross section calculations. Also, the neutron number density 
estimated for this model is somewhat low. 

Schwarzschild and Harm (1967) have recently carried numerical 
calculations through the first four million years of shell helium burning 
for a Population II star of 1.2 solar masses. The thermal instability 
associated with this shell burning phase results in "relaxation oscilla- 
tions" in -which the main helium shell flash leads to the formation of a 
convective zone proceeding outward from the shell. When, after nine 
such cycles, the convective zone reaches the hydrogen-rich layers above, 
the mixing of protons with the products of helium burning can proceed. 

Sanders (1967) has considered the implications of this behavior 
with regard to s-process nucleosynthesis. He finds that for relatively 
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small admixtures of protons into a carbon (oxygen) region, the major 

]_Q ]7 

product will be J C( 0). For both of these cases (a,n) reactions 

proceeding on a somewhat longer time scale may release sufficient neutrons 

to drive the slow neutron capture process. 

These conditions are promising as well with regard to the formation 

of the anomalous abundance peaks discussed above. The mixing of protons 

Ik 

into a carbon and oxygen region will result in the production of N 
as well as and ^0; in fact, for large proton admixtures the 

13 17 

abundance will exceed those of C and 0. In the presence of an 

14 

appreciable helium abundance , this N may be destroyed by 

■ L Sj(o',Y ) i8 F(P + v ) l8 0(a ,Y ) 22 Ne 

The subsequent passage of a shock wave through such a region may result 
in the production of a substantial neutron flux by 

pO pR pP 

Ne (a , n ) Ofe (or , n ) Si 

These neutrons would be produced in a region in which slow neutron 
capture may previously have taken place, providing seed nuclei. We are 
currently undertaking a study of this process. 
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Figure Captions 

Figure 1. Evolution in the density-temperature plane of selected 

representative mass zones expelled by the 2 model of the core 
of a massive star. Details of the calculation are given in 
Arnett (1967a). Radiation pressure dominates above the line so 
designated, and electron degeneracy pressure dominates below the 
line so designated. 

Figure 2. The temperature and density profile employed in our network 
calculations. This density curve corresponds to the low adiabat; 
for the high adiabat, the density is reduced by a factor of ten. 

Figure 3- Schematic of the nuclear reaction network employed in these 
calculations . The reactions involving the light nuclei are 
indicated. Above 12 C, a linear chain through ^ 2 Ni is defined by 
(n,v) and (0,7) reaction links. 

Figure 4. The evolution of a region composed initially of neutrons and 
protons in the ratio N/P = 8 is shown for the low adiabatic cooling 
curve. The mass fractions of the more abundant constituents are 
plotted as a function of time. 

Figure 5- The evolution of a region composed initially of neutrons and 

protons in the ratio N/P = 8 is shown for the high adiabatic cooling 
curve. The mass fractions of the more abundant constituents are 
plotted as a function of time. 
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Figure 6. The evolution of a region composed initially of neutrons and 

protons in the ratio N/P = 3 is shown for the high adiabatic cooling 
curve. The mass fractions of the more abundant constituents are 
plotted as a function of time. 

Figure 7- The evolution of a region composed initially of neutrons and 

protons in the ratio N/P = 8 is shown for the high adiabatic cooling 
3 3 

curve: P = 7x10 T n . The mass fractions of the more abundant 

y 

constituents are plotted as a function of time. 

Figure 8. The equilibrium abundances of nuclei at the freezing temperature 
(T = 4xl0 9 °K) are shown as a function of mass number for the high 
density adiabatic cooling curve. 

Figure 9* The equilibrium abundances of nuclei at a temperature 
9 o 

T = 3x10 K for the high density adiabatic cooling curve are shown 
as a function of mass number. 

Figure 10. The equilibrium abundances of nuclei at the freezing temperature 
(T = 4xl0 9 °K) are shown as a function of mass number for the low 
density adiabatic cooling curve. 

Figure 11. The number densities of neutrons, protons and alpha particles 
and the total number density of heavy elements (A > 4) are plotted 
as a function of temperature for the two adiabats considered. 

Figure 12. Trends among heavy nuclei of abundances attributable to neutron 
capture on fast and slow time scales and to bypassed processes. 
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